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Abstract: We study the Gregory-Laflamme instability of a large uniform black brane 
wrapping a two-sphere compactification manifold. This paper continues the work [1], 
where the compactifications on p-torus were considered. The new features of the spherical 
case are the non-zero curvature of the compactification manifold and the absence of the 
rescaling symmetry due to a built-in stabilization mechanism. We calculate the order of 
the phase transition in dependence on the number d of extended dimensions using the 
Landau- Ginzburg approach. It is found that for d > 11 a uniform spherical black brane in 
microcanonical ensemble exhibits a smooth second order phase transition towards a stable 
branch of non-uniform black brane solutions. The critical number of extended dimensions, 
for which there is a change in the order of the phase transition, is different for microcanon- 
ical and canonical ensembles and does not coincide with the critical number of dimensions 
in the case of the flat toric compactifications. We briefly discuss the origin of this mismatch 
in the orders of phase transition for the different ensembles. 



1 Introduction and summary 



Uniform extended black branes in the presence of compact extra dimensions are unstable 
with respect to long wavelength perturbations if the black brane horizon size is substantially 
smaller than the size of the extra dimensions. This is known as the Gregory-Laflamme (GL) 
instability [2] . Therefore a large black brane should undergo a phase transition once its size 
becomes smaller than a certain critical size. By studying non-uniform perturbations on top 
of the critical black brane it is possible to find the order of this phase transition. It could be 
either a smooth second (or possibly higher) order transition when the black brane becomes 
slightly non-uniform in the compact dimension in a continuos way or a first order transition 
when below the critical size the black brane decays into some completely different solution. 
The final stage of the first order transition is in general unknown (see [3-5] for review). 
The first study of this kind was performed by Gubser for a five-dimensional black string 
on a single compact extra dimension in pure gravity in which case the transition is first 
order [6]. Recently this calculation was generalized by Sorkin to black strings in arbitrary 
number of extended dimensions. He found that the phase transition becomes second order 
in more than twelve extended dimensions [7]. It was also claimed that the critical number 
of dimensions, when there is a change in the phase transition order, depends on whether 
the phase transition happens at fixed black string mass (in microcanonical ensemble) or at 
fixed temperature (in canonical ensemble) . In the latter case the transition becomes of the 
second order in more then eleven extended dimensions [8] . Later Kol and Sorkin considered 
the case with an arbitrary number d of extended dimensions and an arbitrary number p of 
the extra dimensions compactified on the torus T p [1]. In the special case when the sizes 
of all p circles are equal they found that the phase transition order depends only on the 
number d of extended dimensions and not on p. This result is explained by the fact that 
it is thermodynamically preferable for the GL instability to develop only along one of the 
circles on the torus. Therefore the toric black brane with p > 1 behaves effectively like the 
p = 1 black string. 

In this paper we determine the order of the phase transition for a black brane on a two- 
sphere compactification manifold. We use the spontaneously compactified Md x S 2 solution 
of the Einstein-Maxwell theory in D = d + 2 dimensions as the background geometry [9] . 
The presence of the Maxwell field is necessary in order to have a non-flat compactification 
manifold. Aside from having a non-gravitational matter field this case has two important 
features in comparison to the flat compactification set-ups studied earlier which affect the 
properties of the phase transition. First, the two-sphere is not a direct product of two flat 
compact dimensions, and thus it does not support a mode of instability analogous to the 
modes along a single circle on the torus. Instability on S 2 inevitably feels the presence 
of both compact dimensions and in this respect is more similar to the mode on the torus 
when the inhomogeneities along both circles are excited with equal amplitude (this mode 
is referred to as the "diagonal" mode in [1]). The other important difference is that the size 
of the two-sphere is fixed by the parameters of the theory. Therefore there is no rescaling 
freedom which in case of the flat compactifications accounts for the fact that the size of 
the compact dimension can be set arbitrarily. In the terms of dimensional reduction, the 
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radion field is stabilized and has a mass comparable to the inverse radius of the compact 
two-sphere. Because of the absence of an internal length scale for the flat compactification 
case, a set of thermodynamical quantities invariant under rescaling was introduced in [6] 
in order to study the phase transition. The results obtained in such a way correspond to 
the situation when the size of the compact extra dimension is held fixed and the radion is 
infinitely heavy. Since the critical black hole size is comparable to the radius of the two- 
sphere, the radion mass in our set-up is naturally of the same order as the phase transition 
temperature. Therefore one can expect the presence of the dynamical radion to play a 
non-trivial role in the phase transition. The presence of these features suggests that the 
phase transition order for the spherical black brane could be different from the case of the 
flat toric compactification. 

In this paper we follow the method described in [1] to determine the order of the phase 
transition. The method employs the Landau-Ginzburg theory of phase transitions and 
is favoured in comparison to the original Gubser's computation since there is no need to 
compute any third order metric perturbations around the critical black brane. Considering 
the perturbations up to the second order is sufficient in order to compute the free energy and 
the entropy differences between the uniform and non-uniform black brane branches. The 
signs of these differences define the phase transition order in canonical and microcanonical 
ensembles respectively. As a result we find that the transition for the spherical black 
brane in microcanonical ensemble is of the second order when the number of extended 
dimensions exceeds eleven. This coincides with the critical dimension for the "diagonal" 
mode on the two-torus found in [1]. In the canonical ensemble the difference of the free 
energy between the uniform and non-uniform black branes changes its sign when there are 
more than nine extended dimensions. It is lower than the corresponding critical dimension 
for the "diagonal" mode on the two-torus. However, since the canonical ensemble itself 
is ill-defined for the black holes in asymptotically flat extended dimensions due to their 
negative specific heat, it does not make sense to associate this critical behaviour of the 
free energy with the change of the phase transition order. The relation between the order 
of phase transition in the canonical and microcanonical ensembles for a generic system is 
discussed in the appendix A. Therefore we indeed found that the phase transition order for 
the spherical black brane does not coincide with the case of the flat toric compactification. 

The rest of the paper is organized as follows. In section 2 we introduce the set-up and 
briefly describe how to apply the Landau-Ginzburg description of phase transitions to the 
black branes. We provide the details for the perturbative computation of inhomogeneous 
black brane solution in section 3. The results for the phase transition order in various 
number of extended dimension are presented and discussed in section 4. 

2 Set-up and Ginzburg-Landau method 

One of the simplest non-flat compactifications with a stabilized size of the extra dimensions 
is the spontaneous compactification solution of the 6D Einstein-Maxwell theory described 
in a great detail in [9]. In this solution the Maxwell vector field has a magnetic monopole 
configuration on an external two-sphere. It is possible to fine-tune the six-dimensional 
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cosmological constant with respect to the magnetic flux so that the four-dimensional cos- 
mological constant is zero and the remaining four dimensions are flat. It is straightforward 
to generalize this solution to d extended dimensions. The fine-tuning condition between the 
cosmological constant, the vector field coupling, and the Newton's constant in D = d + 2 
dimensions remains the same as in d = 4 case, i.e. 

A ° = -i4r <"> 

The phase transition takes place for a uniform black brane which looks like a spherically 
symmetric Schwarzschild black hole in the d extended dimensions and completely wraps 
the external two-sphere. The Euclidean line element is given by 

ds 2 = f(r)dt 2 + J^dr 2 + r 2 dnj_ 2 + a 2 (d8 2 + sin 2 9 d(f> 2 ) , (2.2) 



where 

/(r) = l-(^)^ 3 . (2.3) 



The first three terms in (2.2) correspond to a Schwarzschild- Tangherlini black hole in d- 
dimensional flat space, and the last two terms correspond to the two-sphere of the compact 
extra dimensions. The radius a of this sphere is determined by the cosmological constant 
scale: 

(2.4) 



,2 



2A D 

The mass and the inverse temperature of such a black brane are related to the horizon 
radius ro by the standard relations for a d-dimensional Schwarzschild black hole 

(d - 2) Q.d-2 ri ^ „ 47r , 

Here Gd = G£>/(47ra 2 ) is the <i-dimensional Newton's constant and £ld-2 is the area of the 
(d — 2)-dimensional unit sphere. 

In the Landau-Ginzburg approach in order to study the order of phase transition in 
canonical ensemble it is sufficient to know the local behaviour of the free energy of the 
system around the critical point. A detailed description of this method in the context 
of black hole phase transitions is given in the reference [1] which we follow closely in our 
calculation. One first computes the local expansion of the free energy as a function of the 
order parameter and the temperature which in our case plays the role of the parameter 
that controls the onset of the transition. The role of the order parameter is played by 
the amplitude A of the inhomogeneous perturbations in the metric and vector field. It 
shows the degree of the deviation of the black brane from the uniform solution and is 
also our perturbative expansion parameter. In order to compute other thermodynamic 
characteristics it is useful to know the free energy as the function of the inverse temperature. 
We expand the free energy up to the fourth order in A around critical point: 1 



F(\;P)~F Q (P) + A[ j-)X z +CX\ (2.6) 



1 We checked that the odd terms in A do not appear in the expansion of the black brane free energy. 
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where Fq((3) = Tq -3 is the free energy for the uniform unperturbed black brane 2 , (3* is the 
inverse critical temperature, and 5(3 = (3 — (3*. The values of the coefficients A and C in the 
expansion (2.6) define the local thermodynamics completely. The phase transition occurs 
when the black brane becomes smaller than a certain critical size, which means that 6/3 
becomes negative. This fixes the sign of A to be positive in order for the uniform phase 
A = to become an unstable extremum of the free energy for 5(3 < 0. In the case when C 
is positive there is a minimum of the free energy for 5(3 < located at 

The presence of the non-trivial minimum in the vicinity of the uniform phase signals a 
smooth second order phase transition towards the slightly non- uniform phase with A = A* . 
The difference in free energies between the non-uniform and uniform black branes is of the 
fourth order in perturbative expansion parameter A and is given by 

W-iM*-£(£) 2 = -c>i. (2-8) 

The mass of the black brane can be obtained from the free energy using the relation 

The entropy S(X; M) can be computed by performing a Legendre transform of (3F{\\ (3) 
with respect to (3 and is given by 

S(A; M) a 5 (M, - ft^ (*f ) »» - ft ( C - 2(< ,_^_ 3) ) A' • (2.10) 

Here Sq(M) = 47rrQ -2 is the uniform black brane entropy, M* is the critical mass, and 
5M = M — M* . The entropy as a function of mass determines the behaviour of the system 
in microcanonical ensemble. In full analogy with the canonical ensemble case the order 
of the phase transition in the microcanonical ensemble is determined by the sign of the 
coefficient in front of the fourth order in A term, which also gives the sign of the difference 
between the entropies of the non-uniform and uniform phases (cf. (2.8)): 

S ^ Sn 1 f C- nfJ , )xi^a 2 Xi. (2.11) 



So d-3 V 2(d-2)(d-3). 
If the non-uniform branch has higher entropy than the uniform one, i.e. when o<i > 0, the 
phase transition in microcanonical ensemble is of the second order, and black brane settles 
in the stable non-uniform branch. 

The free energy of the black brane as a function of the metric and the vector field 
potential is given by the Euclidean action of the Einstein-Maxwell theory evaluated on the 
corresponding solution: 

„ F _ lE ^, vy = (J \ R + [K - A-]) + / fa - A D ) . (2. 12 , 



2 We follow [1] and omit in thermodynamic potentials the constant factor of ^ ^ coming from the GR 



action. 
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Here R is the Ricci scalar, K is the extrinsic curvature on the boundary at infinity and Fn U 
is the vector field strength. The surface term corresponding to a reference geometry with 
the extrinsic curvature K° has to be subtracted in order to make the resulting free energy 
finite. In our case the reference geometry is x M ' -1 x S 2 with the Euclidean time period 
given by the inverse temperature /3 and a fixed radius a of the external two-sphere. The 
configuration space is spanned by the Euclidean solutions for g^ v and that asymptote 
to this reference geometry. 

The non-uniform solution can be found perturbatively by expanding the equations of 
motion and field deviations above the background in the powers of perturbative parameter 
A. The parameter A also plays the role of the order parameter in the free energy expan- 
sion (2.6). In order to compute the free energy up to the fourth order in A it is sufficient 
to find the solution up to the second order. We collectively denote the metric and vector 
field deviations bg^ = g^ - gffl and bV^ = - by X = {bg^ u ,bV^} and expand 
them in powers of A: 

X = A X w + A 2 X (2) + . . . . (2.13) 

The first order perturbation X^ 1 ' is nothing else but the static inhomogeneous Gregory- 
Laflamme mode, and the second order perturbation X^ corresponds to the back-reaction 
of the black brane on the GL mode. 

In practice one expands the free energy in powers of X and then plugs the solution for 
X^ and X^ 2 \ For determining the order of the phase transition one is interested in the 
0(A 4 ) term in the free energy. Using the equations of motion one arrives to the following 
expression for the quartic coefficient C [1]: 

C = F 4 [X^] - F 2 [X^) . (2.14) 

Here the first term is the quartic in the field deviations X term of the free energy expan- 
sion, which is evaluated on the first order solution X^\ The second term i<2[^ 2 )] is the 
quadratic in X term evaluated on the second order perturbation X^ 2 ' . 

3 Perturbative solution 

In this section we are looking for a slightly inhomogeneous black brane solution as a pertur- 
bation around the homogeneous black brane (2.2). The most general ansatz for Euclidean 
metric and vector field, which are static and spherically symmetric in d extended dimen- 
sions reads 

ds 2 = e 2A f(r) dt 2 + ^-dr 2 + e 2C r 2 dU 2 _ 2 + 

fir) 

+ 2a 2 Gdrd9 + a 2 e 2H (e 2J d9 2 + e~ 2J sin 2 9 # 2 ) , (3.1) 

V = — cos6d(j)-aL sin<9# . (3.2) 
2e 

Here the functions {A,B,C,G,H,J,L} = X parametrize the metric and vector field de- 
viations and depend only on the ti-dimensional radial coordinate r and the angles on the 
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external two-sphere and <j). One could, in principle, also include the components of the 
metric, which are odd under the inversion on the two-sphere: 

2 a 2 G dd sin 9 dr d<p + 2 a 2 H dd dO d<j> , 

as well as the even under the inversion component of the vector field (the monopole back- 
ground itself is odd): 

& F even dO . 

However it is always possible to choose a gauge in which these components are not excited 
and, therefore, we omit them. It is useful to expand the set of functions X in spherical 
harmonics on the external two-sphere. It is always possible to choose the direction of the 
inhomogeneous mode to contain only the harmonics with m = that do not depend on <p. 
The spherical harmonics analysis also constrains considerably the possible ^-dependence: 
the first order mode X^ 1 ' contains only the first I = 1 harmonic, whereas the second order 
mode X^ contains only the I = and I = 2 harmonics. Moreover all the functions X can 
be subdivided according to their transformation properties under the coordinate transfor- 
mations on the two-sphere into scalar, vector and tensor ones. The explicit expansion for 
the scalar quantities X s = {A, B, C, H} reads 

X s (r,6) = \xl(r)- cos 6 + \ 2 (x s (r)+x s 2 (r)-^ (3cos 2 #-l)^ , (3.3) 

and the vector quantities X v = {G, L} are expanded as 

X v {r,9) = -Xx\(r) -sin0- X 2 x v 2 (r) -6cos0siii0 . (3.4) 

The tensor harmonic parametrized by J appears only at I = 2 and is given by 

J(r,8) = X 2 h (r) -| sin 2 8. (3.5) 

Hence in this ansatz all the metric and vector field perturbations are parametrized by the 
set of functions xf = {ai,bi,Ci,hi}, x\ = {gi,k}, and j 2 of a single variable r. 

The ansatz (3.1), (3. 2) does not fix the coordinate transformation redundancy in the 
(r, 9) plane and two more gauge fixing conditions have to be specified. These conditions 
can be specified independently at each order of perturbation theory and for each spherical 
harmonic label I. We use this freedom later in order to simplify the resulting equations. 

The non-uniform black brane we are looking for has to satisfy the Euclidean equations 
of motion which in our case read 

R„v = 8tt G d (f^F^ - ^ 9lxv F 2 ^j - 2 - A D , (3.6) 

Vi^ = . (3.7) 

After plugging in the ansatz (3.1), (3. 2) and expanding up to the first order in A one obtains 
a set of ordinary linear differential equations for the functions x\(r) and x\(r). After 
substituting also the background solution conditions (2.1) and (2.4) for the vector coupling 
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and the cosmological constant, the remaining parameters in the resulting equations are the 
radius of the external sphere a and the black brane horizon size Tq. For simplicity we set 
ro to be the unit of length. Then the only parameter in the equations is the dimensionless 
value of the radius a in ro units. A regular solution exists only for a particular value of a 
and defines the static Gregory-Laflamme mode on the black brane. 

In linear order we use the particular choice of the gauge fixing condition proposed 
in [1]: 

_ rfa 1 + 2(d-2)fc 1 
61 " r/' + 2(d-2)/ ' h -°- m 

In this gauge all linear equations are reduced to a single second order equation for c\{r): 



1 U-2 f j\' , 2(d-l)(d-3)/ 



r 



12 



M + jb , 1 , ^2 ci = -» ci . (3.9) 



r d " 2 V V (2(d-2)/ + r/') 2 « 2 

Two conditions have to be imposed on the function c% in order to specify the solution. One 
of them corresponds to the choice of the normalization of c\. We fix it by imposing the 
condition c\(X) = 1 adopted in previous works on black string [1, 6-8]. The other condition 
arises if one requests the solution to be regular at the horizon r = 1 and is given by 

Sr^- (2<i - 2) - (3 - 10) 

However for generic values of the parameter a this regular solution grows exponentially for 
large r. Therefore one is left with a one-parameter shooting problem in which the value of 
a is adjusted so that c\ decays at large r. The equation (3.9) coincides with the first order 
equation (3.6) of [1] for the GL mode in the case of the black string and the black brane in 
the toric compactification. The critical size of the sphere aci at which the GL instability 
sets in is related to the critical wavelength of the GL mode on the black string as: 

k 2 G L = 2/4l • (3-11) 

In the gauge (3.8) all other components of the metric can be explicitly expressed in 
terms of c\{r) as: 



01 = -(d-2)ci , (3.12) 

(ri-2)(r/'-2/) , 1 oN2 f /2^ 

91 = ~ 2r r + W ~ 2) / Cl " r (r f + 2(d - 2) ff {d ~ " ] ^ ~ ^ f + 

+ r(((d - 10)d + 17) / + 2(d - 3)) /' - 4(d - 3)(d - 2)( / - 1) /) c x , (3.13) 
hx = . (3.14) 

At the second order the perturbations contain I = and I = 2 modes which can be 
treated independently. The second order equations are inhomogeneous linear ODE's for 
xq, X2 and 22 that contain source terms quadratic in the first order perturbations found 
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above. We start by considering the zero modes ao,&o,co and ho. The function ho can be 
separated from the other zero modes and is defined by the equation 



'^AoV - ^Aho = Sh^c^gt] . (3.15) 



r d-2 V ' da GL 

The exact form of the source term Sh [bi,ci,gi] is given by (B.l). The regularity condition 
on the horizon fixes the value of the first derivative h' (l) in terms of the value of ho(l), 
while the latter is adjusted in order to find the solution that decays at large r. 

Instead of finding the solution for ao , 60 and cq it is more efficient not to fix the gauge 
at all but to rewrite the equations in terms of the gauge-invariant combinations defined as: 

u = a + 60 - (r co)' , (3.16) 

w = c - ^4 a ■ (3.17) 
r J 

The free energy is obviously a gauge invariant quantity and can be expressed in terms of 
u and w. Thus there is no need to fix any particular gauge and determine ao, 60 and Co- 
The equation for it is a first order differential equation 

u' = ~rh% + S u , (3.18) 

with the source term given in (B.2). This equation can be straightforwardly integrated 
numerically. For the metric to match the flat reference geometry the constant of integration 
u(l) should be chosen so that u vanishes at infinity. The equation for w also happens to 
be of the first order: 

w ' + ^ u+ ^j h °-^ h °-d^h ho = Sw - (3 - 19) 

with the source given in (B.3). The equation can be integrated once the solution for u 
is known. The constant of integration w(l) can be fixed by considering the temperature 
of the perturbed black brane. More precisely the deviation of the temperature from the 
critical one can be expressed in terms of the metric perturbations as 

^ = 6 (1) - ao(l) = u(l) + w(l) + w'(l) . (3.20) 

p* 

Thus different values of w(l) correspond to black branes with different temperatures. How- 
ever, in order to obtain the free energy of the black brane at the critical point we consider 
only solutions with the critical temperature. The integration constant w(l) is then fixed 
by demanding the 5(3 to be zero. The freedom in the choice of w(l) is related to the fact 
that among the linear perturbations around any black hole there always exists a mode 
corresponding to the infinitesimal change of the size of the black hole. 

The equations for the I = 2 spherical harmonic components can also be separated in 
two groups. First we solve for the variables Zi2,j2, and I2, the equations for which are 
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independent from the other variables. We adopt the same gauge condition l 2 = as in the 
linear order. The system of equations for the h 2 and j 2 then takes the form: 



' V- 2 / h' 2 \ = \ C-^h, + 6j 2 )+S h2 , ( 3.2 1 ) 



d -2 y j "i) a 2 \ d 



r 

^=2 (r d - 2 ff 2 )' = -^h 2 +S j2 , (3.22) 

with the source terms given by equations (B.5) and (B.6). It is straightforward to bring 
this system to a solvable from by taking appropriate linear combinations of h 2 and j 2 . 
After that the equations for these linear combinations can be solved separately. In order to 
obtain solution one, as before, has to impose the regularity condition on the horizon and 
solve the one-parameter shooting problem by adjusting the values of the functions on the 
horizon in order to obtain the solution that decays at infinity. 

In analogy to the linear order one can fix the gauge in such a way that it would be 
possible to express the remaining variables a 2 ,b 2 , and g 2 algebraically in terms of c 2 , which 
obeys a single second order ODE. The gauge fixing condition (3.8) is modified at the second 
order by the presence of a term quadratic in the first order variables: 

rfa 2 + 2(d-2)fc 2 4 rf 
62 " rf' + 2(d-2)f " 3r/' + 2(d-2/) h 91 ' (3 ' 23) 

The equation for c 2 is similar to the equation (3.9) that defines the first order perturbation: 

12 



2(d-l)(d-3)f 2 6 



> 1 r d 2 f c'n ) + — — — c 2 - c 2 + 
r d ~ 2 V V ( r // + 2 (d-2)/) 2 a 2 ~ 



(rf' + (d-3)f)f 2 4 
+ 2 (rf + 2{d - 2) f) 2 (k2 + 2J2) + da 2 ^ + ^ = ° ' ^ 

with the source term given by (B.4). Thus finding the solution in this sector boils down 
to another one-parameter shooting problem for the value c 2 (l). After finding the solution 
for c 2 one can determine the remaining variables a 2 and g 2 as: 

, , s (d- l)(d-2) , a 2 r n , 
a 2 = -(d - 2)c 2 -h 2 - 2j 2 - K - g c 2 --fg 2 , (3.25) 

l (d-2)(rf'-2f) 1 1 3rf-2(d-2)f 

92 6 rf' + 2(d-2)f ° 2 12 2 12 rf' + 2(d-2)f 2 
l rf'-2(d-2)f . 

-6rf> + 2(d-2)f j2 + S ^- (3 - 26) 
The explicit expression for the source S g2 can be found in (B.7). 



4 Results and discussion 

Having found the solution for the metric and the vector field we can compute the coeffi- 
cients in the free energy expansion (2.6). By substituting the solution with 5(3 = in the 
Euclidean action (2.12) one obtains an expression for the free energy quartic in A, from 
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d 4 5 6 7 8 9 10 11 12 13 



A 0.544 1.73 3.59 6.15 9.39 13.3 17.9 23.2 29.2 35.8 
C -0.114 -0.384 -0.785 -1.20 -1.41 -1.10 0.14 2.84 7.62 15.2 
a 2 -0.188 -0.317 -0.441 -0.536 -0.576 -0.535 -0.389 -0.112 0.322 0.938 



Table 1. The coefficients A and C in the free energy expansion (2.6) and the entropy variation 
(72 defined in (2.11) for a different number of extended dimensions d. The change of the sign of 02 
between d = 11 and 12 indicates that the phase transition in microcanonical ensemble becomes of 
the second order for d > 11. The analogous change in the free energy behaviour happens between 
d = 9 and 10. 

which the coefficient C can be determined. Alternatively one can use the expression (2.14), 
which is simpler to evaluate numerically. We found the numerical mismatch between the 
two different expressions for C to be less than a percent. By changing the parameters of 
numerical integration we found the change in the values of C to be at the same level, which 
thus can serve as the estimate of the numerical error. The values of C for a various number 
of extended dimensions d are listed in table 1. Note the change of sign of the quartic 
coefficient for d > 9. For a thermodynamically stable system it would mean that the phase 
transition for d > 9 is of the second order in canonical ensemble. 

It is instructive to compare the obtained behaviour of C for a black brane on a two- 
sphere with the case of a flat compactification on the square two-torus T 2 considered by 
Kol and Sorkin in [1]. In the latter case there are two independent inhomogeneous modes 
corresponding to the two circles of T 2 . In order to study the free energy one can consider 
two limiting cases: when only a single mode along one of the two circles is excited, or when 
both of the modes are excited with equal amplitude, the so-called "diagonal" mode [1]. 
The quartic coefficient C in all three cases is presented in dependence of the number of 
extended dimensions d in figure 1. We see that the single direction mode on a torus has 
lower free energy and thus thermodynamically favourable. Due to this fact the toric black 
branes during the phase transition effectively behave like black strings, with only the mode 
along a single circle being excited. In contrast, on the spherical black brane there is only 
one mode, and its dependence on the number of extended dimensions d is different from 
the modes on T 2 . The change of the sign of the coefficient C for the spherical black brane 
happens between d = 9 and 10. 

The behaviour of the black brane in microcanonical ensemble is determined by the sign 
of the coefficient 0*2 in the difference of the entropy between the non-uniform and uniform 
black branes (2.11). In order to find 02 the quadratic coefficient A should be determined. 
There are two independent methods how to compute A. First, one can take a first variation 
of the free energy (2.12) with respect to the temperature, which at the leading order in A 
can be expressed using only the solution with 5(3 = 0. This method was applied in [1] and 
gives 



Alternatively one can use the solution of (3.19) with Sf3 7^ 0, i.e. keep the integration 




(4.1) 
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Figure 1. The quartic coefficient C in the free energy expansion (2.6) for the spherical black 
brane (circles) in a various number of extended dimensions d in comparison to the single direction 
(squares) and "diagonal" (dimonds) modes of the black brane on a square two-torus from [1] . 

constant for the zero harmonic w(l) initially unspecified. The temperature dependence of 
the free energy (2.12) is then obtained by using the relation (3.20) relating 5(3 to w(l). 
The coefficient A is given by 

A = {d-l){d- 2) w(r -> oo) , (4.2) 

where the asymptotic value w(r — > oo) is taken from the solution with 6/3 = 0. 

The resulting values of 02 in dependence on the number of extended dimensions d are 
given in the table 1 and presented in comparison to the case of the toric black brane in 
figure 2. For d > 11 the non-uniform black brane has larger entropy than the uniform one, 
and the phase transition in microcanonical becomes of the second order. We note that 
in microcanonical ensemble the critical number of extended dimensions for the spherical 
black brane coincides with the one for the diagonal inhomogeneous mode of the toric black 
brane. 

For d = 10, 11 the signs of the quartic coefficients in the free energy and entropy, 
and consequently the predicted orders of phase transition in canonical and microcanonical 
ensembles, are different. The possibility of such a situation can be seen already from the 
equation (2.11) for the entropy difference, since the entropy difference can remain negative 
even if the coefficient C would turn to be positive. The details of this effect are discussed 
in appendix A. 
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Figure 2. The coefficient of the entropy difference between the non- uniform and uniform 
spherical black branes (diamonds) in comparison to the single direction (circles) and "diagonal" 
(squares) modes of the toric black brane from [1]. 



A The origin of the phase transition order mismatch in canonical and 
microcanonical ensembles 



We consider a generic form of the free energy expansion in the vicinity of the critical 
temperature ft 



5J3 
ft 



ft 



\P* 



(A.l) 



where A is the order parameter, and 5f3 = (3 — ft, cf. (2.6). The entropy is given by a 
Legendre transform of f3F with respect to f3 and reads 



S(M, A) ~ ft./i + ft<5M + ft 



1 



4(/i + h 



-5M 2 - -, ^ A „ J A4X 2 - ft C- AfJ T, i J A 4 , 

(A.2) 

with SM = M — M* . The specific heat c u of the uniform black brane with A = kept fixed 
is given by 

c u = ~P (^FJ a _ q - - W/i + h) ■ (A.3) 



2(/i + / 2 ) 



A 2 



4(A + / 2 ) 



\PF) \ 



By using this expression the entropy expansion can be rewritten in the following form: 

P*A 2 



S(M, A) ~ ft/i + PJM - ^5M 2 - ^5MX 2 - ft 



A 2 

2c u 



C + 



2c u 



A 



(A.4) 



If the specific heat c u is negative, the coefficient C = C + ^jf- — , which determines the 
stability of the non-uniform branch in microcanonical ensemble, may become negative 
even for positive values of C. 
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In the case of a positive C a non-trivial minimum of the free energy, located at 



appears for 5/3 < 0. Therefore, in canonical ensemble the system should settle in the 
non- uniform phase with A = A*. The specific heat in this non-uniform phase is given by 

(A . 6) 

Comparison with (A. 4) gives the following relationship between the quartic coefficients in 
the two ensembles and the specific heats in the different phases: 



% = — • (A-7) 



This condition holds for any thermodynamic system and tells that the signs of the quartic 
coefficients in the free energy and the entropy expansions are different if and only if the 
specific heat has different sign in different phases. Hence, the mismatch between the phase 
transition orders in canonical and microcanonical ensembles can happen only if the specific 
heat changes its sign in the course of transition from the uniform to the non-uniform 
phase. Such a situation is indeed observed in some gravitational systems (c.f. [12] and 
references therein). If the system is thermodynamically stable before the phase transition, 
i.e c u > 0, then a situation is possible when the non-uniform branch is stable in the 
microcanonical ensemble (C > 0) and unstable in the canonical ensemble (C < 0). In such 
a case the system exhibits a second order phase transition in the microcanonical ensemble. 
In canonical ensemble the transition is of the first order and proceeds towards some third 
phase which is thermodynamically stable. This behaviour is attributed to the fact that in 
this case the system has negative specific heat in the non-uniform phase, as can be seen 
from (A. 7). 

The opposite situation takes place for the spherical black brane in the cases d = 10, 11. 
Specific heat flips its sign from negative on the uniform branch to positive on the non- 
uniform. Nevertheless it does not mean that the branch of the non-uniform branes is 
thermodynamically stable. Due to the change of sign of the specific heat the "stable" 
non-uniform phase, which appears in canonical ensemble for 5/3 < 0, corresponds to the 
black branes with the mass larger than M* . Thus the mass of the brane does not cross the 
GL critical value in the course of this transition, and the would be new phase is related 
to the local minimum of the entropy which appears for 5M > 0. Thus it seems that posi- 
tive specific heat of non- uniform branch in canonical ensemble is spurious, and the phase 
transition in canonical ensemble never proceeds towards the found non-uniform branch. 
Moreover in the case at hand the change of sign of the specific heat during the transition 
would mean that slightly non-uniform spherical black branes with flat asymptotics are 
thermodynamically stable in certain number of extended dimensions which does not seem 
to be the case. 
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B The source terms for the second order perturbation equations 

For the sake of completeness we present here the full expressions for the sources in back- 
reaction equations. The source terms for the I = mode equations (3.15), (3.18) and (3.19): 



q d + 1 /2 2r 1 2 2 r / 1 / d _ 2 ,y 1 

<j« = r Ci ; — o of -oi — j—jj r JCi o< 

3 1 3 fda 2 1 3 / V r ^ ' « 

2d-3 2 d- 1 / /2 2/1 a 2 \ 2 2 / 
^ " -37277 c i " — j7 c i " 3 [jji + J-) 9i + 3 (d _ 2)/ /5i &i+ 

2 1 f r f + (d-3)f 1 \ 2 l r/' + 2(d-2)/ , 

The source term for the C2 equation (3.24) reads 

S C2 = 1 * ( r 2 f' 3 (5a 2 (d -2){d- 1) /' + 4(d - 3)r) + 

2 3a 2 (2(d - 2) / + r f) 3 V V ' 

+ 2(d - 2) / 2 /' (a 2 (d - 6)(d - 3)(d - 1) /' - 4(d(5d - 13) + 3)r) + 

+ (d - 2)r / / /2 (a 2 (d - l)(7d - 27) /' - 2(9d + 8)r) - 8(d - 2) 2 (2d(2d - 5) + 5) / 3 ) c 2 + 



+ 



2(d-2)/((d-3)/ + r/') , 2 2a 2 f 2 ((d-2)f + rf) , 2 
6(d-2)/ + 3r/' Cl + 3(d-2)(2(d-2)/ + r/') 5l + 



2 / (2(d - 6)r //' + 4(d - 2)(d - 1) f 2 + 5r 2 f 2 ) , 

H 9 c i <7i + 

3(2(d-2)/ + r/') 2 



+ - , 2a 2 dr 2 /' 4 + 2(d - 2)r / /' 2 (4o 2 d /' + r) + 

3(d-2)d(2(d-2)/ + r/') V V 7 

+ (d - 2) / 2 /' (a 2 d(7d - 15) /' - 4(d + 4)r) - 16(d - 2) 2 (d + 1) / 3 ) g 2 + 

, 2/((15d-32)r//' + 2(d-2)(8d-15)/ 2 + 5r 2 /' 2 ) , 

H c i 3i+ 

3(2(d-2)/ + r/') 2 

+ 3r(2( ,_ 2 ) / + r/f 2 (^ " ^ >* + ^ " 2 ^ " ^ / /' 2+ 
+ (d - 2)(7(d - 5)d + 40)r / 2 /' + 2(d - 5)(d - 2) 3 / 3 ) c[ c x + 

2/((3d-5)r/' + (d-5)(d-2)/) „ , 2a 2 f 2 ((d - 2) / + r /') „ 



+ 



6(d - 2) / + 3r /' Cl Cl + 3(d - 2) (2(d - 2) / + r /') 9l 91 



2a 2 / (6(d - 2)r f f + (d - 2) 2 / 2 + 4r 2 f 2 ) , 

3(d-2)r(2(d-2)/ + r/') + 



+ ^2 (3(d-4)r 3 // 3 + 2(d- l)(4d- 7)r 2 f f' 2 + 

3r(2(d-2)/ + r/') 3 ^ 

+ 2(d - 2)(d(9d - 26) + ll)r f f + 8(d - 2) 3 (d - 1) f + 3r 4 / /4 ) ci ffl . (B.4) 
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The source terms for the rest of I = 2 mode equations (3.21), (3.22) and (3.26): 

S h = J±^*> c\-lA~\f^\^ (r^f 9l )' bl + lf 9l * , (B.6) 

- - (d "^" 2) r/ , + 2 ( d _ 2)/ (2/^^-2/(^ + 3/4^) " 

1 f 7r ^ 2 r r , 7rf + 8(d-2)f h 
~ r/' + 2(d-2)/ (l8? 6l + 18 75151 + 3 /Ml + 18 h 91 



1 a 



2 



36 rf + 2{d-2)f 



rf(9 2 i)" + {9rf + 4(d-2)/) / 5l <7i+ 



+ 1 (5, /' + 2(d - 2) /) /' rf) - 1 r// + ; 7 _ 2)/ ■ (B.7) 
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